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ABSTRACT 

Context. Studying the physical conditions in circumstellar disks is a crucial step toward understanding planet formation and disk evo- 
lution. Of particular interest is the case of HD 100546, a Herbig Be star that presents a gap within the first 13 AU of its protoplanetary 
disk, a gap that may originate in the dynamical interactions of a forming planet with its hosting disk. 

Aims. We seek a more detailed understanding of the structure of the circumstellar environment of HD 100546 and refine our previous 
disk model that is composed of a tenuous inner disk, a gap and a massive outer disk (see Benisty et al. 2010b). We also investigate 
whether planetary formation processes can explain the complex density structure observed in the disk. 

Methods. We gathered a large amount of new interferometric data using the AMBER / VLTI instrument in the H- and K-bands to 
spatially resolve the warm inner disk and constrain its structure. Then, combining these measurements with photometric observations, 
we analyze the circumstellar environment of HD 100546 in the light of a passive disk model based on 3D Monte-Carlo radiative 
transfer Finally, we use hydrodynamical simulations of gap formation by planets to predict the radial surface density profile of the 
disk and test the hypothesis of ongoing planet formation. 

Results. The SED (spectral energy distribution) from the UV to the millimeter range, and the NIR (near-infrared) interferometric data 
are adequately reproduced by our model. We show that the H- and K-band emissions are coming mostly from the inner edge of the 
internal dust disk, located near 0.24 AU from the star, i.e., at the dust sublimation radius in our model. At such a short distance, the 
survival of hot (silicate) dust requires the presence of micron- sized grains, heated at ~ 1750K. We directly measure an inclination 
of 33° + 11° and a position angle of 140° + 16° for the inner disk. This is similar to the values found for the outer disk (/ ^ 42°, 
PA ^ 145°), suggesting that both disks may be coplanar. We finally show that 1 to 8 Jupiter mass planets located at ~ 8 AU from the 
star would have enough time to create the gap and the required surface density jump of three orders of magnitude between the inner 
and outer disk. However, no information on the amount of matter left in the gap is available, which precludes us from setting precise 
limits on the planet mass, for now. 

Key words. Stars: individual: HD 100546 - Stars: circumstellar matter - Techniques: interferometric 



1. Introduction 

Studying the physical conditions in circumstellar disks and the 
processes that rule the evolution of gas and dust provides the 
context for the formation of planets. It is now clear that sev- 
eral processes are simultaneously acting on the disks at a given 
time. With the advent of powerful mid- and far-infrared space 
telescopes such as Spitzer, a new class of objects has been iden- 
tified, the p re-transitional and transitional disks (e.g., |Espaillat| 
et al.|2010| ). Pre-transitional disks have a typical spectral energy 
distribution (SED) with a near-infrared excess resulting from the 
emission of hot dust and gas located in an inner disk, a dip in the 
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* Based on observations collected at the VLT (ESO Paranal, Chile) 
with programs 082.D-0010(A), 083.C-0298(A,B), 083.D-0224(C), 
083.C-0146(B), 083.C-0144(A,D), 083.C-0236(A), 075.C-0637(A) 



mid-infrared range likely caused by a gap, and at longer wave- 
lengths, the signature of an optically thick outer disk. 

HD 100546, well studied already by ISO, is a late Herbig Be 
star (B9.5Ve) that is now classified as a transitional disk. Located 
at ~103^g pc (as measured by Hipparcos, [van den Ancker et alT] 
,1998 ), it has a large-scale disk inclined by ~4U" (U" meanmg 
pole-on, this convention will be used throughout the paper) ex- 



tending up to 4" that was first imaged in scattered light ([Pantin 



et al., 



2000 ). Its estimated age of 10 Myr ( van den Ancker et al. 



1997| ) corresponds to t he timescale o n which disks are found 
to dissipate ( Hillenbrand et al.||2QQ8| ), which suggests that the 
disk is evolved. Coronagraphic imaging revealed a large-scale 
envelope (-1000 AU) and a disk extending u p to 515 AU with 
an asymmetric brightness profile ( Augereau et al .||2001[ Grady 



|et al.||200T] ). Based on its SED that shows a dominant mid 
infrared excess, peaking at 40 pm and with a weak 2-8 pm emis- 
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sion, Bouwman et al. (2003 ) suggested the presence of an inner 
hole within 10 AU that could possibly result from the presence 
of a very low-mass stellar or planetary c ompanion in the disk 
(20 Mjup, Acke & van den Ancker|2006 ). Nulling interferome- 
try at mid-mtrared wavelengths showed that the radial tempera- 
ture law in the few inner AUs cannot be explained with a con- 
tinuously flared disk ( |Liu et aL 2003 1 ). Additional evidence of 



a cavity is provided by means of spectroscopy: HST/STIS ob- 
servations of H2 eni ission revealed a hole extending to 13 AU 
( |Grady et all2005at , w hile a lack of 6300 A [01] emi ssion in- 
side 6.5 AU is observed ( |Acke & van den Ancker|2006 ). Recent 
CO ro-v ibrational observations ( Brittain et al.|2009 van der Plas| 



et al. 2009 ) reported an extended emission with a large inner hole 
inside -11-13 AU, confirming the status of pre -transitional disk 
given to HD 100546. 

Recent far-infrared and millimetric observations enabled fur- 



ther st udies of the gas and dust in the outer disk. Panic et aL] 
( |2010| ) detected pure rotational lines of CO and derived a mass 
of 10~^ Mq of molecular gas in the disk. Their model is con- 
sistent with an inclined outer disk (i~40°) in Keplerian rotation 
that extends out to 400-500 AU in radius. The Herschel Space 
Observatory detected several additonal molecular lines as well 
as forbidden atomic lines ([CII] and [01]) that probably resu lt 
from photo-dissociation by stellar photons ( [Sturm et al.|2010| ). 

On the other hand, very little is known about the inner 
disk. Far ultraviolet (FUV) observations of molecular hydro- 
gen ([Lecavelier des Etangs et al. [2003 [ Martin-Zai'di et al!] 



2008|) imply that the observed gas lies very close to the cen- 



tral star. Although these observations give evidence that the in- 
ner 1.5 AUs are not cleared of gas, clues are still missing for a 
firm conclusion about the location of this gas (FUV-driven pho- 
toevaporative wind from the surface layers of the disk, gas in 
an inner disk...). However, the near-infrared (NIR) excess, the 
emission from gas lo cated as close to the central star as 0.5 AU 
(Acke & van den Ancker 2006), as well as the signs of accretion 
( Deleuil et al.|2004] rall suggest that the innermost AUs are not 
empty. 

This was confirmed by interferometric measurements in the 
K-ba nd that spa tially resolved the dust emission at the sub-AU 
scale ( [Benisty et al. 2010, hereafter BIO). The authors, based on 
two visibility measurements at 2.2 and 8.7 yum, proposed a disk 
model that includes an inner disk, a gap, and an optically thick 
outer disk that are able to reproduce the SED and the visibilities. 
See their Fi g. 3 for a sketch of the disk configuration. Benisty 
2010| considered a stellar luminosity of 22 L© to repro- 



et al. 



duce the photometric measurements at the shortest wavelengths, 
assuming a null Ay and a normal 7?v=3.1 extinction law, while 
other authors previously used a higher value of 32 L© using a 
non-zero Ay but an anomalous extinction law, 7?v~5.1 ( |van den| 
[Ancker et al.lT997] ). 

Owing to the small number of interferometric measurements 
used in BIO, some questions remained unsolved. What is the 
exact structure of the inner disk? Is it asymmetric and puff'ed 
up, or rather smooth? Is there any hot matter located inside the 
dusty rim that could emit at shorter wavelengths? Are the inner 
and outer disks coplanar? Could the presence of a sub- stellar 
companion be constrained? 

In this paper, we present new spectro-interferometric mea- 
surements obtained in the K-band, as well as the first spatially 
resolved measurements in the H-band and first closure phase 
measurements for HD 100546. We build a complete SED from 
previously published photometric and spectroscopic data, and 
reassess the stellar luminosity in our model. Based on the ex- 
tensive dataset, we revisit the disk morphology and study the 



eff'ects of a putative planetary companion on the disk structure 
and evolution. In § 2 we present the observations and the data 
processing. In § 3 we model the SED and interferomeric data. 
In § 4 we use hydrodynamical simulations to predict the surface 
density profile in the disk when a planet carves a gap into it. § 5 
is the summary. 



2. Observations and data processing 

2.1. Interferometry 

HD 100546 was o bserved at th e Very Large Telescope 
Interferometer (VLTI; [Sch511er|2007[ ), using the AMBER instru- 
ment that allows the simultaneous combination of three beams in 
the near-infrared (Petr ov et al.[[2007| . It delivers spectrally dis- 
persed interferometric observables: visibilities, closure phases, 
and diff'erential phases. In the following, we present H- and K- 
band observations taken in the low spectral resolution mode (LR; 
R~30) with the 1.8 m Auxiliary Telescopes (ATs). The data 
were obtained during several Guaranteed Time Observation pro- 
grammes. HD 100546 was observed with 11 diff'erent baselines 
of four VLTI configurations during eight nights from January to 
May 2009. The longest (projected) baselines achieved for the H- 
and K-band observations are ~92 and -127 m, respectively, cor- 
responding to a maximum angular resolution of A/B ~3.9 mas 
and 3.6 mas, respectively. We use the VLTI nomenclature to 
identify the diff'erent configurations. A summary of the ob- 
servations used in this paper is given in Table [T] In addition 
to HD 100546, calibrators (HD104479, HD100901, HD77049) 
were observed to correct for instrumental eff'ects. Their proper- 
ties can be found in Table [2] All observations were performed 
without the fringe-tracker FINITO. 

The data reduc tion was performed f ollo wing standard proce - 
dures described inl Tatulli et aL] ( [2007] ) and [Chelli et"aLl ( [20091 ), 
using the amdlib package, release 3.0b4, and the yprick inter- 
face provided by the Jean-Marie Mariotti Centei[^ Raw spec- 
tral visibilities and closure phases were extracted for all the 
frames of each observing file. A selection of 80% of the high- 
est quality frames was made to overcome the eff'ects of instru- 
mental jitter and unsatisfactory light injection. Calibration of the 
AMBER-hVLTI instrumental transfer function was done using 
measurements of the calibrators, after correcting for their diam- 
eter. Finally, calibrated visibilities and closure phases were av- 
eraged over the wavelength range of each band (H: from 1.5 to 
1.8 yum; K: from 2.0 to 2.4 yum), to obtain broad-band observ- 
ables. Because of the combined eff'ects of an intrinsically lower 
flux of the source and of a lower instrumental sensitivity, the per- 
formance of AMBER is worse in the H-band than in the K-band. 
As a consequence, our H-band data are of lower signal-to-noise 
ratio, which translates into two efl'ects: (i) the visibility points 
show a larger scatter in the H-band, and (ii) we were unable 
to estimate the H-band closure phase for the longest baselines 
where the source is almost fully resolved. The lower signal-to- 
noise ratio of the H-band data also explains the lower number 
of data points in the lower panel of Fig. [T] especially at long 
baseline lengths where the visibilities are lower. 

2.2. Spectroscopy 

FUSE spectrum: Although HD 100546 was observed t wice with 
FUSE, in 2000 and in 2002 (Martin-Zaidi et al.'2008), we only 
considered the exposures of 2002 that have the best signal-to- 
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Table 1. Log of the observations and associated (u, v) coverage. All data were taken in 2009. 
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Table 2. Stellar and calibrator properties. 



Star 


K 


H 


Sp. Type 


Diameter [mas] 


HD100546 


5.4 


5.9 


B9Vne 




HD101531 


4.2 


4.3 


Kim 


0.6+0.1 


HD104479 


3.9 


3.9 


KOIII 


0.8+0.1 


HD100901 


4.1 


4.2 


KO/KIIII 


0.8+0.1 


HD77049 


4.4 


4.6 


K2III 


0.7+0.1 



noise ratio. The observations were all made in the time-tagged 
mode, using the 3C x 3C LWRS aperture at a resolution of 
about 15 000. The total exposure time was split into several sub- 
exposures that were co-aligned using a linear cross -correlation 
procedure and added segment by segment. The individual spec- 
tra (905 - 1187 A) were processed with the CALFUSE pipeline 



2000 Dixon et al.l 



processing software v3.0.7 (Sahnow et al. 
|2Q07| ), to correct for instrumental effects. 

UVES spectrum: The star was observed on 2005 March 21, 



with UVES ( [Dekker et al.|2Q00^ ), a cross-dispersed echelle spec- 
trograph mounted at the VLT. For these observations, the blue 
arm was used with a spectral coverage between -3700 A and 
5000 A. We used no image slicer, low gain with a 1x1 binning of 
the CCDs and a slit width of 0.4 providing the maximum resolu- 
tion of -80,000. Wavelength calibration frames were taken with 
a long slit and a ThAr arc lamp. The spectrum w as reduced using 
the UVES pipeline v3.2.0 (Bal lester et al.|2Q0Q| ) available on the 
ESO Common Pipeline Library. The spectrum was corrected for 
the Earth's rotation to shift it in the heliocentric rest frame. No 
standard star is available to provide an accurate flux calibration. 
We scaled the UVES flux to the lUE one, which agrees well with 
the best photospheric model for HD 100546 (see Sect. |3.1| ). 

lUE spectra: We used archived lUE spectra at short (SWP: 
1150 - 1975 A) and long (LWP: 1850 - 3350 A) wavelengths to 
cover the entire UV band. The spectra, of high signal-to-noise 
ratio, were obtained with the large aperture (10'' x 20'') and 
the high-dispersion mode which produces a two-dimensional 
echelle spectrum containing approximately 60 orders, with a res- 
olution of roughly 0.2 A. 

2.3. The spectral energy distribution 

The photometric data used in the SED are the same as in 
BIO. They were retrieved from the literature. They include pho- 



tometry from ground-based telescopes in the optical and NIR, 
2MASS d ata, mid-infrared photometry from the IRAS and ISO 
satellites ("Ardila et al. 112007; Henning et aLl|1994t |Grady et al." 



2001 ; Malfait et al. 1998; Sturm et al. 2010])). These are shown 
as black dots in Fig. [2] We also included spectra from the FUSE, 
lUE, and ISO archives. 



3. Modeling 

3.1. Stellar properties 

Various estimates of the fundamental parameters of HD 100546 
have been proposed over the past 15 years. The efl'ective tem- 
perature found in the literature varies from 10000 to 11 000 K, 
estimates of the visual extinction toward this star range from Ay= 
0.0 to Av= 1.03, and the stellar luminosity estimates range from 



22 to 32 Lq (e.g., |l 


rhe etal.|1994l 


van den Ancker et al.|1998l 


Valenti et al.||2000t 


Acke & van den Ancker||2006K These sig- 


nificant diff'erences 


lave motivated a thorough modeling of the 



photosphere. We started considering photospheric models from 



the Kurucz grid of LTE models ([Kurucz 1979; Castelli & Kurucz 



2004 ). The sub-grid thus obtained samples the parameter space 
9000 K < Teff < 12000 K with 250 K steps and 3.5 < log ^ < 
5.0 with 0.25 dex steps. Then, syntheti c spectra were calculated 
in NLTE using the Synspec package of |Hubeny & Lanz ( 1995 ). 
This code allows for an approximate NLTE treatment of lines 
with LTE photospheric model^ The result of this stage is to de- 
termine the LTE model that provides the best fit to the observed 
spectra. The model assumes a solar composition, a helium abun- 
dance, Y = He/H = 0. 1 by number, and a microturbulent velocity, 
= 2 kms"^ For details of the modeling method of FUV and 
UV spectra of Herbig stars, we refer the reader to Bo uret et al.| 
( 2003) . The best-fit model provides the following stellar proper- 
ties: Teff = 10500 K, log g = 4.0 and v sin / = 55 km s"^ . 

The stellar luminosity was calculated from T^^ and using 
HIPPARCOS parallaxes ( van den Ancker et al. 1998; Bertout| 
et al. 1999|) and bolometric c orrections using the method de- 
scribed in |BesseU et al.| ( [T998] ): 



log 



5 X log(l/7i) -my- BCy - 0.26 
Z5 



(1) 



where n is the parallax in arcseconds, my is the V-band magni- 
tude of the object, and BCy, the bolometric correction. We found 
a luminosity of 25.9^2 9 ^o- previously published values are 



http://nova.astro.umd.edu/Synspec43/ 



3 



Tatulli et al.: The inner disk of HD 100546 



K-band 




id and position angle PA^. The ring thickness is fixed at 20% of 
the radius. The modeled visibility can be written as 



Effective Baseline (m) 

Fig. 1. Broad-band squared visibilities with (Icr) error bars 
shown against eff'ective baseline for the K-band (upper panel) 
and the H-band (lower panel). The H-band plot shows less data 
points because of a lower number of measurements above the 
required S/N ratio threshold. The predictions of the best uniform 
ring model are overplotted in solid line. The filled blue points 
correspond to the visibilities used in the BIO paper. 



within 2cr of our estimate. In the following, we use 26 . We 
note that BIO used a stellar luminosity of 22 L© and a null 
extinction, assuming a normal (Ry= 3.1) extinction law. In ad- 
dition, they had a limited sampling of the stellar flux at wave- 
lengths shorter than 3000 A. A comparison with the stellar mod- 
els of'S iess et al.| ( |200 ) suggests that this luminosity is probably 
too low, unless the metallicity of HD 100546 is Z=0.01 (half so- 
lar). 

The dataset used by BIO can also be fitted properly with a 
higher luminosity, e.g., 32 L© as in van den Ancker et al. (1997]) 
or with 26 L© (this work) but with a non-zero and anomalous 
extinction (Ay- 0.2 - 0.3 and Ry- 5.5). The anomalous extinc- 
tion law is needed to fit the flat UV spectrum now included in 
the SED. We note that the photospheric model and extinction 
curve we use lead to a slight underestimation of the UV flux near 
lOOnm. We suggest this excess is a consequence of the accretion 
process still going-on for HD 100546. This has no significant 
impact on the continuum energy budget and continuum emis- 
sion from the dust disk because the flux densities involved are 
too low by a large factor. 



V 



[H,K] 



ring 



T/ Z7[^'^] _i_ T/ zr[H,K] 
j^[H,K] j^[H,K] ' 



(2) 



with K = 1 (unresolved central star) and Flf^^/(Fl"^^^+Fl^^^^) 
the NIR excess, i.e., the excess over the stellar flux divided by the 
total flux in the H- and K-band respectively. Modelling the stellar 
photosphere emission as described in Sect. |3.1[ we estimate H- 
and K-band excesses of ^ 0.53 and ^ 0.75. 

An independent fit of the H- and K-band visibilities for 
HD 100546 shows that the continuum emission in both bands 
comes from the same location, within the uncertainties. We 
therefore performed a classical x'^ minimization using the com- 
plete H- and K-band visibility datasets together. We constrain the 
NIR-emitting region to be located at Rd = 0.25 AU ± 0.02 AU, 
which reinforces the estimation made in our previous study 
(Rd = 0.26 AU; BIO). We note that this estimation is indepen- 
dent of the continuum emission mechanism, it is in particular 
independent of the dust properties in the disk. 

Because the observations provide a good coverage of the 
spatial frequencies (see Table [T]) we are also able to estimate 
the inclination and position angle of the inner disk. They are 
id = 33° ± ir and PAd = 140° ± 16° respectively. The best 
ring model is shown in Fig.[T] together with the squared visibil- 
ities plotted against the eff'ective baseline length. The eff'ective 
baseline length is defined as 



^eff = B ^|cos^{e) -h cos^(id)sin^(6). 



(3) 



where 6 is the angle between the baseline direction and the major 
axis o f the disk, and id is the disk inclination ( [Tannirkulam et al.| 
2008 ). This representation is useful to show a large dataset in a 
concise way once the inclination and position angles of the disk 
are known. Our estimates, relevant for the inner disk, agree well 
with the previously published values of both the inclination and 
position angles of the disk measured at larger scales (e.g., |Ardila 
|et al. 2007] ). This result suggests that the inner and outer disks of 
HD 100546 are most likely coplanar - as implicitly assumed in 
B 10 - or very close to it. This is unlike the case of HD 135344 B 
for example, another transitional disk (Grady et al.^2009j . In the 
rest of this paper we will use a single inclination and position 
angle for the whole disk, adopting the values commonly used in 
previous modeling, namely 42° and 145°, respectively. 

We remark that if our fit is good for the long (^ 30m) eff'ec- 
tive baselines (i.e., for angular resolution better than ~ 12 mas), 
our simple ring model is systematically too high with respect 
to the visibilities at short baselines. Although most of the NIR 
emission is located around 0.25 AU, this may indicate that there 
is some extra extended emission not included in the ring model. 
The origin of this extended emission will be discussed in the 
next section, together with our radiative transfer disk model that 
describes the disk environment of HD 100546. 



3.2. Location of the near-infrared emission 

To derive the basic characteristics of the NIR emission, we first 
consider a simple geometrical model made of a ring of uniform 
surface brightness to roughly estimate the location of the emis- 
sion. We choose a ring instead of, e.g., a uniform disk because 
the bulk of the emission is expected to come from the inner rim 
of the dusty disk, as shown in BIO. The ring model has three 
degrees of freedom, namely the ring inner radius Rd, inclination 



3.3. Disk structure 

Radiative transfer with MCFOST: Going beyond the simple ge- 
ometrical ring model, we used the Monte Carlo-based 3-D radia- 
tive transfer code MCFOST (Pinte et al.,2006„2009 ) to compute 
the SED and NIR images of HD 100546. Visibilities and closure 
phases are computed from the NIR model images. Here, we re- 
fine the model developed in B 10 that was based on one visibility 
measurement. In particular, we use a higher stellar luminosity of 
26 Lo (see, Sect.[37T]). 
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100.0 



Wavelength (|Linn) 

Fig. 2. Observed spectral energy distribution of HD 100546 
compared to our best-model prediction (full black line, 
Teff=10500 K, 26 L©, Av=0.2, Ry=5.5). The photosphere is 
shown as a red line. The photometric data used in B 10 are repre- 
sented by black dots. The binned FUSE spectrum is plotted with 
red dots, the binned lUE spectrum with blue dots, the archival 
ISO data are plotted with cyan dots for ISO-PHOT, green dots 
for ISO-SWS, and magenta dots for ISO-LWS. The IRAS-LRS 
data are represented by yellow dots. 

Table 3. Best-model parameters. Boldface indicates when the 
value differs from the BIO model. 
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1 
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Note: /3 is the exponent of the scale height radial profile, p is the expo- 
nent of the surface density radial profile. The grain size distribution is a 
power-law with a slope of -3.5. 



The disk is similarly composed of (i) an inner dusty disk with 
an inner edge at 0.24 AU, (ii) a gap between 4 AU and 13 AU, 
necessary to reproduce the dip in the SED, (iii) a surface layer 
(over the outer disk) of small grains, and (iv) an outer disk from 
13 AU to 500 AU holding most of the mass (see the sketch of 
the disk in Fig. 3 of BIO). It is difficult to provide reliable error 
bars for all these values without a full exploration of the param- 
eter space. But the uncertainties, or validity ranges, in particular 
regarding the position and size of the gap, were discussed in B 10 
(see §5). Their discussion is still valid. 

Considering our larger dataset and the change in the stel- 
lar luminosity, some of the BIO model parameters need to be 
slightly adjusted. The location of the inner rim is slightly de- 
creased from 0.26 to 0.24 AU to better match the visibilities. 
We also reduced the mass of both the inner disk and of the sur- 
face layer (see Table [3]) to reproduce the NIR excess as well 
as the IMIR bump in the SED. The outer disk now extends up 
to 500 AU, i nstead of 350 AU, following recent CO line ob- 
servations by Panic et al. ( 2010| ). We modified the outer disk 
mass accordingly. In order to generate the large silicate emission 



feature seen in the SED at lO ^um, we used amorphous olivine 
grains (Dorsc hner et al.p995 l for the composition of the sur- 
face layer, whe reas the rest of the d isk is made of astronomical 
silicate grains ( [Draine & Lee|p^84 ), as in BIO. The remaining 
parameters such as density profile, scale height, flaring and grain 
sizes were kept unchanged. The parameters are summarized in 
Table [3] This model provides an adequate fit of the SED from 
the UV to the millimeter range, as shown in Fig.|2] 

Fig.|3]compares the model predictions with the observations. 
The models agree well with the entire set of H- and K-band vis- 
ibilities and closure phases. The "kink" in the model visibility 
curves at B~10m is caused by the outer disk that has a sharp 
inner edge and with a finite sampling of the baselines in the 
models. Increasing the sampling makes the curve smoother, but 
the "kink" remains. A further demonstration of this can be seen 
in Fig. 4 (top row, middle panel), where the "kink" disappears 
when only the inner disk is considered in the field-of-view of the 
interferometer. 

In addition to the visibilities that probe the inner rim location 
and the properties of the dust it contains, we now discuss the first 
closure phase measurements of HD 100546. The closure phase 
is related to the level of asymmetry of the source, here of the NIR 
continuum emission. Roughly, the values of the closure phases, 
given in radians, give the ratio between the asymmetric and sym- 
metric fluxes for the corresponding telescope triplet (Monnier 
r2000 ). As theoretically expected (La chaume 2003 }, the depar- 
ture from centro- symmetry becomes significant for the closure 
phase for the longest baselines only, i.e., when the source is well 
resolved. In our case, with B^ax = 127 m we find a closure phase 
signal of ~ 5 - 10° in the K-band, whereas it is compatible 
with 0° in the H-Band. Recall, however, that the uncertainties 
are larger and there is a lack of visibilities at very long baselines 
atH. 

Properties of the inner disk: From the data we constrain the 
location of the disk inner edge at r ~ 0.24 AU from the star. 
Assuming silicate grains, the inner disk is composed of micron- 
size particles whose albedos make the scattered light an impor- 
tant part (~ 30%) of the disk total K-band emission, as dis- 
cussed in BIO. Fairly large grains need to be used (i.e., with 
<^max = 5 yum) to make them survive (i.e., not sublimate too fast) 
at the inner edge of 0.24 AU, a value imposed by the interfer- 
ometric observations. Because of the higher stellar luminosity 
than the one used in BIO, we find a temperature for the hottest 
grains of ~ 1750 K. We note that this value is on the warm side 
for silicate grains but we did not attempt to use diflferent compo- 
sitions. We did not explore either the possiblity that much larger 
grains be present in the inner disk because their contribution to 
the disk emission in the NIR is not dominant with respect to the 
grains used here. Because of the small NIR excess in the SED, 
we find that the inner disk must be tenuous, with a dust mass of 
^ 1.75x10-10 Mq. 

We find an integrated radial optical depth of r ^ 10 (in the 
I-band) in the equatorial plane of the inner disk. Because it is 
opaque, depending on its geometrical thickness it may (or may 
not) cast a shadow and may (or not) significantly modify the 
temperature structure of the outer disk, hence its scale height. 
We verified that the scale height profiles used in our parametric 
models agree with a disk in hydrostatic equilibrium. To check 
this, we calculated the hydrostatic scale height of a vertically 
isothermal disk, consistent with the vertical Gaussian density 
profiles we used, and used the two extreme temperatures found 
at each radii to bound the scale height of the actual disk. As a 
minimum temperature for each radius we used the temperature 
of the disk midplane, the maximum temperature being estimated 
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Fig. 3. K-band (top) and H-band (bottom) MCFOST modeling of HD 100546. From left to right: NIR images (normalized intensity 
to one at maximum, in logarithmic scale), visibility (red solid lines) and closure phase (red crosses), compared with the interfero- 
metric observations (black circles and error bars). In the middle panels, the "kink" in the model visibility curves at B~10 m is a real 
feature caused by the sharp inner edge of the outer disk. 



at the disk surface. We find that the scale height of the outer disk 
(i.e., H=12 AU at r=100 AU) falls between the two estimates, 
which suggests that the inner disk does not dramatically modify 
the structure of the outer disk. We note that the scale height of 
the inner disk (6 AU at 100 AU) is slightly too large, by about 
30%. 

The measured closure phase value of ~ 5 - 10° in the K- 
band is well reproduced by our model. It is a direct consequence 
of the anisotropic scattering of starlight by the dust in the inner 
disk. The phase function of the dust grains in our model has an 
asymmetry factor g ~ 0.6 (g = is isotropic scattering and 
^ = 1 is fully forward-scattering). The forward- throwing nature 
of the light scattering process in our model directly results in the 
asymmetry of the NIR brightness distribution because the disk 
is inclined. Note also that the closure phase signal is entirely 
caused by the inner disk, because the outer disk is too faint to 
contribute at K-band, as shown in Fig. |4] (middle). The outer disk 
has a diff'erent impact on the visibilities, see below. We remark 
here that we do not need to invoke the presence of a "pufFed-u p" 
inner rim (e.g., 'Dullemo nd et al.|2001t|lsella & Natta|2005| ) to 
reproduce the interferometric data because the scale height of 
the inner rim follows the smooth law defined for the whole inner 
disk, i.e., H{r) = Ho(r/ro)^. In our model we thus have H(Ri^) = 
0.015 AU corresponding to a ratio H/R of 0.06 for the inner rim, 
this is a factor 2 to 4 smaller than the values of H/R~ 0.1 - 
0.25 that are typically expected when the rim is puff'ed-up and 
vertical ( [Dullemond et al.|20QT] ). |Monnier et al.| ( |2006| ) showed 
that sharp "puff'ed-up" inner rims are expected to produce large 



closure phase signals as soon as they are resolved enough. In 
our case, the value of the inner rim diameter (0.48 AU) in units 
of fringe spacing A/Bmax is ~ 1.3, a value for w hich substantial 
closure phases of ^ 40° would be expected (see [Monnier et al. 
,2006^ Fig. 22). Such high values are clearly ruled out by our 
data. A sharp "puff'ed-up" inner rim for HD 100456 is unlikely. 

The closure phase signal is also very sensitive to the value 
of the flux ratio between the star and the disk. Because the ra- 
tio will vary with the stellar luminosity, a careful determination 
of this parameter is essential. The study of the closure phase is 
a go od way to cross-check our luminosity estimation of Sect. 
3.1 Figure [4] (left) shows the behavior of the (visibility and) clo- 
sure phase for = 22 L© , 26 L©, and 32 L©, including the 
values used in the literature. While the visibility profile does 
not change significantly with the luminosity, the closure phase 
increases from a few degrees to 15° as the stellar luminosity 
decreases. Such analysis indirectly confirms that the luminosity 
of HD 100546 must indeed He in the [22,32] L© range. More 
accurate data would be of great interest to reduce this range of 
validity. 

Gap and outer disk: The Auxiliary Telescope's field of 
view of A/D ^ 250 mas = 25 AU enables us to probe the outer 
environment of HD 100546, that is, the gap and the outer disk. 
This extended emission within the field of view directly impacts 
on the visibility profile at short baselines, as shown in Fig. |4] 
(middle and right). Including the emission from the outer disk in 
the model and the field of view better reproduces the visibilities 
at short-baselines - despite their scatter - than the ring model 
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Fig. 4. K-band MCFOST modeling of visibility (top) and closure phase in degrees (bottom) for various values of model parameters. 
Left: Various stellar luminosities are compared: 22 L© (green), 26 L© (red), and 32 L© (blue). Middle: Effect of the field-of-view 
when using the AT (red) or UT (black) telescopes. With the latter, only the inner disk of HD 100546 is seen, which is shown 
in green. The black and green curves lie on top of each other, as expected. Right: the interferometric observables are plotted for 
different locations of the inner edge of the outer disk, 2 = 9 AU (blue), 13 AU (red), and 17 AU (green). 



alone. Typically, with short baselines of a few tens of meters, 
the interferometer becomes mostly sensitive to the location of 
the inner edge of the outer disk, as emphasized by Fig. [4] (right). 
This suggests that the gap must end somewhere between 9 AU 
and 17 AU. Unfortunately, the scatter in the data does not allow 
for a better estimation of this parameter. 



Interestingly enough, the outer part of the disk would not 
have been traced with the 8 m Unit Telescopes of the VLTI be- 
cause it lies outside its field-of-view (A/D ^ 60 mas = 6 AU). 
This points out the importance of the auxiliary telescopes, not 
only in terms of adding baseline capabilities to that of the UTs, 
but also in terms of field-of-view for tracing potential extended 
emission. Indeed, lower visibilities at short baselines with re- 
spect to the unifor m ring case were als o measured for other 
Herbig AeBe stars ( Monnier et al.|[2006] ) and were interpreted 
as resulting from an extended halo. Our model shows that the 
scattered light emission from the outer part of a flared disk (even 
when continuous), a contribution that is usually neglected when 
modeling the K-band emission, can naturally and easily explain 
this fairly common feature. 



4. Hydrodynamical simulations of gap formation 
and inner disl< evolution 

It is well documented that a planet in a disk is expected to carve 
a gap, a trench, in its density profile. Numerous studies exist 
to study the form ation of these gaps by planets via numerical 
simulations (|Bryden et al.[T999||Kley|1999| [Nelson et al.|2000| 



Varniere et al.||2004| ), but none to our knowledge has focused 
on trying to reproduce the density jump values needed to match 
the observations. In this paper, our goal is not to revisit the gap 
formation mechanism but rather to verify whether a planet em- 
bedded in the disk of HD 100546 can explain the large jump 
in the surface density profile observed on both sides of the gap. 
To do so, we carried out a series of 2D simulations involving 
planets of various masses and opening disk gaps. We have made 
several assumptions and simplifications: circular orbit for the 
planet, constant accretion rate onto the star over the its full life- 
time and equal to the present day accretion rate, the accretion 
rate in the disk is the same as onto the star and the same at all 
radii in the disk. However, we have run the numerical simula- 
tions sufficiently long to l et the inner disk evolve until ste ady- 
state is reached. Following Ac ke & van den Ancker (2006), the 
lowest planet mass estimate based on the local disk thickness is 
about 1 Jupiter mass, we therefore started our study with a one 
Jupiter-mass planet. 
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Fig. 5. Left: evolution of the density as function of radius in the case of a 6 Jupiter-mass planet. The solid line indicates the density 
profile used in our model. Right: final density profile as a function of the disk radius for various planet masses. 



To reproduce the density jumps, several mechanisms are 
competing in the inner region. These are (i) The accretion onto 
the star; for HD 100546 we have an estimated accretion rate of 
~ 10"^ Mo/yr ([ Peleuil et al.|200 4HGra dy et al.|2005b). This will 
tend to deplete the inner region, which in turn needs to be re- 
plenished. Indeed, with the estimated mass available in the inner 
region and the observed accretion rate, the inner region should 
empty itself in less than one year, (ii) The planet "pushes" the 
gas in the inner region during the creation of its gap. This is a 
transitory event and occurs soon after the planet formation, it is 
not part of the steady state we are looking for. (iii) Some material 
from the outer disk can flow accross the gap, cross the planet's 
orbit, and replenish the inner region. 

To reproduce the density jump suggested by the observations 
of HD 100546, one needs to reproduce a surface density for the 
inner disk that is about three orders of magnitude lower than the 
surface density just outside of the gap. Ideally, it also needs to 
produce the proper gap width, i.e., between about 4 and 13 AU. 



4.1. Hydrodynamic code and simulation setup 

To cover a sufficiently long time- span in the simulation we used 
the algorithm FARGO ( |Masset| [2^0} |2002 |), which efiminates 
the azimuthal velocity from the computation of the Courant- 
Friedrich-Levy condition. This speeds up the computation and 
facilitates the study of the long-term disk evolution. Then, to 
follow the disk over several orders of magnitude in radius to see 
how the density evolves close to the st ar as well as outside, we 
used the code FARGO- 2D ID ( Crida et al. 2007), which is an ex- 
tension of the standard version of FARGO, where the 2D grid is 
surrounded by a simplified ID grid made of elementary rings 
that are non azimuthally resolved. It aims at studying and taking 
into account the general evolution of the disk in the simulations 
of planet-disk interactions. 

We ran several simulations to find and optimize the planet 
mass needed to create the required density jump. A planet 
will open a gap near its Lindblad resonance located at = 

e) where m is the mode of the Lindblad res- 



onance, rp the position of the planet and e the eccentricity of the 
orbit. To obtain the outer edge at ~ 13 AU, we placed the planet 
at 8 AU. We set up the 2D disk around it to capture the 2D struc- 
ture of the wave. The 2D disk runs from 4 to 20 AU and the ID 
disk extends from 0.1 AU to 150 Alj^to reproduce the global 
behavior of the whole region. 

4.2. Simulated density profile: Evolution with time 

In Fig. [5] (left) we follow the time evolution of the surface density 
in the first 100 AU of the disk for a 6 Jupiter-mass planet. Such a 
planet produces a surface density jump of the correct amplitude: 
the inner disk surface density drops by three orders of magnitude 
in about 200 000 years, confirming that a 6 M/^^^ planet located 
at 8 AU can carve the disk of HD 100546 with the proper surface 
density jump across the gap. 

Figure |5] (right) shows that most planet masses are able to re- 
produce the strong density jumps imposed by the observations. 
The only differences are the time required to reach that state, the 
amount of mass left in the gap, and the width of the gap. The 
main parameters that regulate how much mass is left in the inner 
region are the accretion rate onto the star and the disk viscosity. 
If we can estimate the accretion rate from observation, it is a 
lot harder to have access to the local disk viscosity. We thus 
adopted an ad hoc value of the viscosity that was successfully 
used for previous simulations (see for example Quil len et al.| 
^2005 ) and aim only at reproducing the density jump. This does 
not allow us to provide a lower limit for the mass of the putative 
planet because the longest timescale is well within the age of 
the system and so far we do not have constraints on the mass 
inside the gap. 

Further remarks can be made regarding these numerical 
results. It is interesting to note that the density profile of the 
inner region becomes fiat early on (before 10 000 years) owing 



^ We ran test simulations with different sizes of both disks to vali- 
date the range of the 2D vs ID part of the domain and found no major 
differences. 
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to the accretion onto the star. It becomes flatter than the r"^ 
power-law used in our model. However, as mentioned above, 
the slope of the density profile is not well constrained by the 
current observations. We verified that radiative transfer models 
that use an inner disk with a flat inner density profiles (i.e., 
constant with radius) also match the interferometric data and 
SED. 

Also, during the 700 000 years of the simulation time span 
the planet has a slow migration from 8 to about 7 AU. This 
slightly moves the outer edge of the gap but it remains within 
the observational limit. However, in all simulations the gap width 
appears too narrow compared to the surface density profile esti- 
mated from the data. We did not try to solve this question. We 
can only suggest that an elliptical orbit for the planet may help 
to broaden the gap. Also, HD 100546 is significantly older than 
the final time of each simulations. It is therefore difficult to es- 
timate how much wider the gap would have grown in the simu- 
lations by the age of the star. Finally, we also note that the full 
extent of the inner disk, in particular its outer radius, are not 
well contrained by the data. Rout of the inner disk, which defines 
the inner boundary of the gap, could be moved outward some- 
what depending on the exact vertical density profile. We refer the 
reader to a forthcoming paper for a more detailed comparison of 
the profiles from the hydrodynamical simulations and those from 
data fitting. For now, we only note that a planet can produce the 
observed surface density jump across the gap. 

5. Summary 

The main results of our paper are summarized below: 

- using the AMBER/VLTI interferometer in H- and K-band, 
we spatially resolved the near-infrared emission region of 
HD 100546. Most of this emission arises from the inner edge 
of its inner disk located at 0.24 + 0.02 AU from the star 

- a puff'ed-up rim in the inner disk appears unlikely because of 
the small closure phase signal 

- for the first time, our observations also constrain the in- 
clination and the position angle of its inner disk. We find 
/ = 33° ± 11° and PA = 140° ± 16°, which suggests that the 
inner and outer disks are likely coplanar. 

- we described the circumstellar environment of HD 100546 
in the light of a passive disk model based on 3D Monte- 
Carlo radiative transfer. Our model is composed of (i) a ten- 
uous inner disk with a dust mass of ~ 1.75x10"^^ M© from 
0.24 AU to ~4 AU, (ii) a gap devoid of dust, and (iii) a mas- 
sive outer disk with a dust mass of ~4.3xl0~^ M© extending 
from -13 AU to -500 AU. The model reproduces both the 
SED and the interferometric observations adequately. 

- we show from hydrodynamical simulations that disk-planet 
interactions with a planet located at ~8 AU are able to repro- 
duce the observed density jump between the inner and outer 
disks, across the gap. 
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